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The {P + l3 + ) Qv and (e/? + ) 0l , modes of 96 Ru, 102 Pd, 106 Cd, 124 Xe, 130 Ba and 156 Dy isotopes are 
studied in the Projected Hartree-Fock-Bogoliubov framework for the + — > + transition. The relia- 
bility of the intrinsic wave functions required to study these decay modes has been established in our 
earlier works by obtaining an overall agreement between the theoretically calculated spectroscopic 
properties, namely yrast spectra, reduced B(E2:0 + — » 2 + ) transition probabilities, quadrupole mo- 
ments Q(2 + ) and gyromagnetic factors g(2 + ) and the available experimental data in the parent and 
daugther even-even nuclei. In the present work, the required nuclear transition matrix elements are 
calculated in the Majorana neutrino mass mechanism using the same set of intrinsic wave functions 
as used to study the two neutrino positron double-/? decay modes. Limits on effective light neutrino 
mass (m„) and effective heavy neutrino mass (Mjv) are extracted from the observed limits on half- 
lives Ty 2 (0 + — > + ) of (/3 + /3 + ) „ an d { £ P + ) 0l/ modes. We also investigate the effect of quadrupolar 
correlations vis-a-vis deformation on NTMEs required to study the (/3 + /3 + ) 0i/ and (e/3 + ) 0t/ modes. 

PACS numbers: 23.40.Bw, 23.40.Hc, 27.60.+J, 27.70.+q 



I. INTRODUCTION 

The sixteen rare, experimentally distinguishable, 
modes of nuclear f3(3 decay, namely the double-electron 
emission (/3 - /3 - ), double-positron emission ([3 + f3 + ), 
electron-positron conversion and double-electron 

capture (ss) with the emission of two neutrinos, no neu- 
trinos, single Majoron and double Majorons, arc scmilcp- 
tonic weak transitions involving strangeness conserving 
charged currents. The (3 + (3 + , e(3 + and ee modes are en- 
ergetically competing and we shall refer to them as e + (3(3 
decay. The experimental as well as theoretical study of 
nuclear (3~ [3~ mode has been excellently reviewed over 
the past decades, which can be found in the recent review 
[l[ and references there in. Also, the experimental and 
theoretical studies devoted to the e + BB decay have been 
reviewed over the past years @, H, 13) S B 0, II, H, 03 EH • 
Owing to the confirmation of flavour oscillation of neu- 
trinos at atmospheric, solar, reactor and accelerator neu- 
trino sources, it has been established that neutrinos have 
mass. However, it is generally agreed that the obser- 
vation of {/3/3) 0v decay can clarify a number of issues 
regarding the nature of neutrinos, namely the origin of 
neutrino mass (Dirac vs. Majorana), the absolute scale 
on neutrino mass, the type of hierarchy and CP viola- 
tion in the leptonic sector, etc. Further, the possible 
mechanisms for the occurrence of the lepton number vi- 
olating ((3(3) 0u decay are the exchange of light as well as 
heavy neutrinos and the right handed currents in the 
LRSM, the exchange of sleptons, neutralinos, squarks 
and gluinos in the i? p -violating MSSM, the exchange of 
leptoquarks, existence of heavy sterile neutrinos, compos- 
iteness and extradimensional scenarios. In nine Majoron 
models, namely IB, IC, IIB, IIC, IIF, ID, IE, IID 



and HE [Tj], the single Majoron accompanied neutri- 
noless double beta {f3(3cj)) 0l/ decay and double Majoron 
accompanied neutrinoless double beta ((3(3(jxf)) 0u decay 
occur in the former five and the latter four, respectively. 
The study of {(3(3) 0u decay can provide stringent limits on 
the associated gauge theoretical parameters and its ob- 
servation can only ascertain the role of various possible 
mechanisms in different gauge theoretical models. 

In principle, the f3~ f3~ decay and e + (3(3 decay can 
provide us with the same but complementary informa- 
tion. The observation of (e + (3/3) 2u decay modes will be 
interesting from the nuclear structure point of view, as 
it is a challenging task to calculate the nuclear transi- 
tion matrix elements (NTMEs) of these modes along with 
((3~ /3~),„ mode in the same theoretical framework. Fur- 
ther, the observation of (e + (3[3) Qv decay modes will be 
helpful in deciding issues like dominance of mass mech- 
anism or right handed currents [l3j |. In an attempt to 
study the role of m v , A and rj mechanisms, Klapdor- 
Kleingrothaus et al. have analyzed the 71.7 kg.y data 
collected from 1990-2003 on enriched 76 Ge [3 and have 
shown that there is an apparent degeneracy in the param- 
eters [lf|. It has been also concluded that the analysis 
of a high sensitive (f3~ I3~) 0v experiment e.g. 76 Ge and 
a suitable high sensitive mixed mode decay e.g. 124 Xe is 
more advantageous 

In spite of the fact that the kinetic energy release in 
the (se) Qv mode is the largest, the experimental and the- 
oretical study of this mode has not been attempted so 
far. The conservation of energy-momentum requires the 
emission of an additional particle in the {ee) Qv mode. 
Further, the emission of one real photon is forbidden for 
the + — * + transition if atomic electrons are absorbed 
from the iv-shell. Therefore, one has to consider various 
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processes such as internal pair production, internal con- 
version, emission of two photons, L-capture etc. [(|. The 
decay rates of the above mentioned processes have to be 
calculated at least by the third order perturbation the- 
ory. Resultingly, there is a suppression factor of the order 
of 1CP 4 in comparison to the {e(3 + ) 0u mode. Hence, the 
experimental as well as theoretical study of (e + (3(3) Ql , de- 
cay has been restricted to ((3 + (3 + ) Ql/ and (g/3 + ) n ,, modes 
only. Arguably, Sujkowski and Wycech [16[ have shown 
that there will be resonant enhancement of the (se) Qv 
mode if the nuclear levels in parent and daughter nuclei 
are almost degenerate i.e. Q — (E 2P — E 2 $) ~ 1 keV, 
where the energy difference is for atomic levels. Inter- 
estingly, Barabash et al. have reported that there might 
be a degeneracy between the 112 Sn ground state and an 
excited 0+ state at 1870.9 keV in 112 Cd fulfilling the res- 
onance enhancement condition for the {se) 0v mode [l7| . 
It is expected that the study of this (ee) 0l/ mode may be 
interesting in the near future. 

The complex structure of nuclei in general, and of mass 
region 96 < A < 156 in particular, is due to the subtle 
interplay of pairing and multipolar correlations present 
in the effective two-body interaction. The mass regions 
A ~ 100 and 150 offer nice examples of shape transi- 
tions at N = 60 and 90, respectively. The nuclei are 
soft vibrators for neutron number N < 60 and N < 90 
and quasi-rotors for N > 60 and N > 90. Nuclei with 
neutron numbers N — 60 and 90 are transitional nuclei. 
The yrast spectra of Te and Xe isotopes, on the other 
hand, follow an approximate inverse parabolic type of 
systematics with minimum energy of 2 + states occurring 
for 120 Te and 120 Xe isotopes, respectively. In this mass 
region 96 < A < 156, the deformation parameters (3 2 are 
in the range (0.1409 ± 0.0046) - (0.3378 ± 0.0018) corre- 
sponding to 132 Xe and 156 Gd isotopes, respectively and 
hence, it is clear that deformation plays a crucial role 
in reproducing the properties of these nuclei. In nuclear 
(3(3 decay, the role of deformation degrees of freedom in 
addition to pairing correlation has been already stressed 
Recently the effects of pairing and quadrupolar 
correlations on the NTMEs of {(3~ /? _ ) „ mode has been 
studied in the ISM @,[S|. In the PHFB modeb the role 
of deformation effects due to quadrupolar [22l 123 . [U, [25| 
and multipolar correlations [26| has been also studied. 

The shell model is the best choice for calculating the 
NTMEs as it attempts to solve the nuclear many-body 
problem as exactly as possible. However, the first ex- 
planation about the observed suppression of M 2v was 
provided in the QRPA model by Vogel and Zirnbauer 
|27T | and Civitarese et al. [28| . Further, the QRPA and 
its extensions have emerged as the most successful mod- 
els in correlating single-/? GT strengths and half-lives of 
((3~ [3~) 2u mode. In spite of the spectacular success of 
the QRPA in the study of (3(3 decay, the necessity to in- 
clude the deformation degrees of freedom in its formalism 
led to the development of the deformed QRPA model for 
studying (3(3 decay of spherical as well as deformed nu- 
clei. The effect of deformation on the [J3~ (3~) 2v mode 



for the ground state transition 76 Ge — > 76 Se was stud- 
ied in the framework of deformed QRPA with separable 
GT residual interaction [29j and, very recently, employ- 
ing realistic forces [3(| ■ A deformed QRPA formalism to 
describe simultaneously the energy distributions of the 
single-/? GT strength and the ((3~(3~) 2l/ mode matrix ele- 
ments for 48 Ca, 76 Ge, 82 Se, 96 Zr, 100 Mo, 116 Cd, 128 < 130 Te, 
136 Xe and 150 Nd isotopes using deformed Woods-Saxon 
potential and deformed Skyrme Hartree-Fock mean field 
was developed |3l| . Rodin and Faessler [32| have stud- 
ied the (3~(3~ decay of 76 Ge, 100 Mo and 130 Te isotopes 
and it has been reported that the effect of continuum on 
the NTMEs of {(3~ (3~) 2lJ mode is negligible whereas the 
NTMEs of ((3~(3~~) 0u mode are regularly suppressed. 

In the PHFB model, the interplay of pairing and de- 
formation degrees of freedom are treated simultaneously 
and on equal footing. However, the structure of the in- 
termediate odd Z-odd N nuclei, which provide infor- 
mation on the single-/? decay rates and the distribu- 
tion of GT strengths, can not be studied in the present 
version of the PHFB model. In spite of this limita- 
tion, the PHFB model, in conjunction with pairing plus 
quadrupole-quadrupole {PQQ) [33j interaction has been 
successfully applied to study the + — > + transition of 
{(3~ (3~) 2v mode, where it was possible to describe the 
lowest excited states of the parent and daughter nuclei 
along with their electromagnetic transition strengths, as 
well as to reproduce their measured (3~ (3~ decay rates 
[12, [24]]. The main purpose of using the PQQ interaction 
is to study the interplay between sphericity and defor- 
mation. In this way, the PHFB formalism, employed in 
conjunction with the PQQ interaction, is a convenient 
choice to examine the explicit role of deformation on 
the NTMEs. The existence of an inverse correlation be- 
tween the quadrupole deformation and the size of NTME 
M 2v has been also confirmed [12, [H, [24| ■ In addition, it 
has been observed that the NTMEs for [3~{3~ decay are 
usually large in the absence of quadrupolar correlations. 
With the inclusion of the quadrupolar correlations, the 
NTMEs are almost constant for small admixture of the 
QQ interaction and suppressed substantially in realistic 
situation. It was also shown that the NTMEs of (3~ (3~ 
decay have a well defined maximum when the deforma- 
tion of parent and daughter nuclei are similar and they 
are suppressed for a difference in deformations in agree- 
ment with previous QRPA calculations [29| . The defor- 
mation effects are also of equal importance in the case of 
{(3-p-) 2v and [P~(3-) 0v modes {MM- 

Moreover, the PHFB model along with the PQQ in- 
teraction in conjunction with the summation method has 
been successfully applied to study the (e + (3(3) 2 decay 
of 96 Ru, 102 Pd, 106 ' 108 Cd, 124 > 126 Xe, 130 ' 132 Ba [l| 
and 156 Dy [34[ isotopes for the + — > + transition, 
not in isolation but together with other observed nuclear 
spectroscopic properties, namely yrast spectra, reduced 
B(E2:0 + — > 2+) transition probabilities, quadrupole mo- 
ments Q(2 + ) and gyromagnetic factors g(2 + ). This suc- 
cess of the PHFB model has prompted us to apply the 
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same to study the + — > + transition of (/3 + /3 + ) 0l/ and 
(e/3 + )o„ modes for the above mentioned nuclei. It has 
been observed that in general, there exists an anticorre- 
lation between the magnitude of the quadrupolar defor- 
mation and the NTMEs M 2v of (e + f3(3) 2v decay. In the 
case of (e + /3/3) 2l/ decay, we observed that the deforma- 
tion plays an important role in the suppression of M^v by 
a factor of 2-13.6 approximately [HI, [24], H3|- Therefore, 
we aim to study the variation of NTMEs of (/3 + /3 + ) „ 
and (£/3 + ) „ modes vis-a-vis the change in deformation 
by changing the strength of the QQ interaction. 

The present paper is organized as follows. The theoret- 
ical formalism for calculating the half-lives of (f3 + f3 + )n l , 
and (e/3 + ) 0l/ modes has been given by Doi et al. Q. 
Hence, we briefly outline steps of the detailed deriva- 
tions in Sec. II. In Sec. Ill, we present the results and 
discuss them vis-a-vis the existing calculations done in 
other nuclear models. In the study of (/?/?) „ decay, the 
practice is to either extract limits on various gauge theo- 
retical parameters from the observed limits on half-lives 
of the (/3/?) „ decay or predict half-lives assuming cer- 
tain value for the neutrino mass. Presently, the avail- 
able experimental limits on half-lives of (/3 + /3 + ) „ and 
(e/3 + ) 0t , modes arc not large enough to provide strin- 
gent limits on the effective gauge theoretical parameters 
(m„) and (Mm)- Therefore, we also predict half-lives 

T l/2( 0+ -> 0+ ) 0f {P + & + )ou aild ( £ P + )ov m0deS for 96Ru > 

102 Pd, 106 Cd, 124 Xe, 130 Baand 156 Dy isotopes, which will 
be helpful in the future experimental studies of {e + Pf3) 0v 
decay. In addition, we study the deformation effect on 
NTMEs of {P + (3+) 0u and (e/3 + ) „ modes and show that 
the NTMEs have well defined maximum for similar de- 
formations of parent and daughter nuclei and they are 
suppressed for a difference in deformations. Finally, the 
conclusions are given in Sec. IV. 



II. THEORETICAL FORMALISM 

In the Majorana neutrino mass mechanism, the effec- 
tive charged current weak interaction Hamiltonian den- 
sity Hw for (3 + decay due to VF-boson exchange including 
hadronic currents can be written as 



Hw — 



G_ 



The left handed V — A leptonic and hadronic currents 
for P + decay are given by 



v eL^i ii (1 - 75) e, 



31 = 



(2) 
(3) 



where g v — cos9 c and 8 C is the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing angle for the left and right 
handed d and s quarks. Further, 
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U ei N lL . 
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The Majorana neutrino field JVj has mass rrii and the 
mixing matrices U of left handed neutrinos are normal- 
ized i.e. Wei\ 2 = 1- 

i 

Usually, the decay rates for the + — > + transition of 
(/3 + /3 + ) „ and (£(3 + ) 0l/ modes are derived by making the 
following assumptions: 

(i) The light and heavy neutrino species of mass mi < 10 
eV and rrii > 1 GeV, respectively are only considered. 

(ii) The nonrelativistic impulse approximation is as- 
sumed for the hadronic currents. 

(iii) The recoil current is neglected. However, it has been 
shown by Simkovic et al. [35] and Vergados [36j that 
the consideration of pseudoscalar and weak magnetism 
terms of recoil current reduce the NTMEs up to 30%, 
which needs to be further investigated. 

(iv) The S1/2 waves describe the final leptonic states. 

(v) The calculation of phase space factors is made easier 
by considering no finite de Broglie wave length correction. 

(vi) The CP conservation is assumed. Consequently, the 
effective light neutrino mass (m„) and effective heavy 
neutrino mass (Afjv) are real. 

With these approximations, the inverse half-lives Ty 2 

for the 0+ -> 0+ transition of (/3 + /5 + ) 0iy and (e/3 + ) 0l , 
modes in 2n mechanism are given by @ 



^1/2 00) 



G01 {P) (M GT - Mp) + 



G01 09) {M GTh - M Fh y 



m e J \ (M N ) J 



(M N ) 

gt - M F ) (M G Th - M F h) , 



05) 



where /? denotes the (/3 + /3 + ) „ / { £ P + )q u mode and 

(m v ) = V' U%rm, rrii < 10 eV, (6) 

£ 

(Mn)- 1 = Y"u 2 el mr\ mi > l GeV . (7) 



In the closure approximation, NTMEs Mp, Mqt, Mph 
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and Mqtii are written as 



Mp 



f E(°fII^k-t-||o+) 



= bn-V m H(r)T-T-\\0+), 



(8) 
(9) 



M Fh 
M GTh 



MMpme)- 1 (£?-) ^(0+||5(r)r-r-||0+) 
4tt (M p m e ) _1 (o£ ||cr„ • cr m <5 (r) t~t~ || Oj) . 

i 



(10) 

(ii) 



The neutrino potential H (r) arising due to the exchange 
of light neutrino is defined as 



H(r) 



4nR 
(2^ 



d 3 q 



exp (iq • r) 
w (w + ~A) ' 



with 



A=(E 



{Ej + E F ) 



(12) 



(13) 



In addition, the inclusion of effects due to finite size 
of nucleons (FNS) and short range correlations (SRC) is 
required. The FNS is usually taken into account by a 
dipole type of form factor making the replacement 



9v -> gv 



A 2 



A 2 + fc 2 



and 



gA ->■ gA 



A 2 



A 2 



k 2 
(14) 

with A = 850 MeV. In the PHFB model, the con- 
figuration mixing takes care of the long range correla- 
tions. The effect of short range correlations (SRC), which 
arise mainly from the repulsive nucleon-nucleon poten- 
tial due to the exchange of p and lj mesons, is usually 
absent. To study the (/3~/3 - ) 0l , mode, the SRC has 
been incorporated by Hirsch et al. through the ex- 
change of w-meson 1371 , Kortelainen et al. [38| as well 
as Simkovic et al. [39j by using the unitary correlation 
operator method (UCOM) and Simkovic et al. [40j by 
self-consistent CCM. This SRC effect can also be incor- 
porated through phenomenological Jastrow type of cor- 
relation using Miller and Spencer parametrization by the 
prescription 



3?3?J\fOf\XXJ), (15) 



where 



f(r) 



'(1-br 2 ) (16) 



with a = 1.1 fm~ 2 and b = 0.68 fm~ 2 [il[. It has been 
shown by Wu and co-workers [42| that for the (/3 _ /? _ ) „ 
mode of 48 Ca, the phenomenologically determined /(r) 
has strong two nucleon correlations in comparison to the 
effective transition operator fOf derived using Reid and 
Paris potentials. 

In the PHFB model, the calculation of the NTMEs 
M a (a = F, GT, Fh and GTh) of the (/3+/3+) 0v and 
(e/3+) 0iy modes is carried out as follows. The two basic in- 
gredients of the PHFB model are the existence of an inde- 
pendent quasiparticle mean field solution and the projec- 
tion technique. To start with, amplitudes (ui m ,Vi m ) and 
expansion coefficients Cy ;m required to specify the axi- 
ally symmetric HFB intrinsic state |<&o) with K — are 
obtained by carrying out the HFB calculation through 
the minimization of the expectation value of the effec- 
tive Hamiltonian. Subsequently, states with good angu- 
lar momentum J are obtained from |<i>o) using the stan- 
dard projection technique [43| given by 



8tt 2 



D^ (n)R(n)\^ )dn, (17) 



where R(fl) and D^ Q (U) are the rotation operator and 
the rotation matrix, respectively. Further, 



\$o)=Y[(u lm + v tm bl m bl rfi )\0) 



(18) 



with the creation operators b\ and i>| defined as 



t>im — '^^^ia.mO^am a nd b m — / ](— m Ci a . m a} a _ r 
a a 

(19) 

Finally the NTMEs M a of the (/3 + /? + ) „ and (e/3+) „ 
modes are given by 
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M a = (^-°||O a r-r-||* 



J t =0\ 
00 / 



,Ji=0Jf=0 1-1/2 



/ ' n (Z M),(z-2,N+2){0) ^2 (af3\O a T r | jS) 
o 

E 



I Au)* \ 

\J Z-2,N+2)eP 



(u)* 

2,N+2 



1 + ^Z,N^)fz-2M+2 



■ sin #<i#, 



(20) 



where 



{det[l + F«(0)/ (,r)t ]} 1/2 
x{det[l + J , M(e)/M t ]} 1 / a dJ (ff) sin(0)d0 



(21) 



and 



n(z,A0,( Z - W a)(«9 - {det[l + F^(0)/Wt !jv+2 ] } i/2 

xldetll + F^^)/^.^]} 1 ^. 

(22) 

The 7r(^) represents the proton (neutron) of nuclei in- 
volved in the (/3 + /3 + ) 0ly / (e/3 + ) 0l/ mode. The matrices 
/z,Af and Fz,n{6) are given by 

(23) 



and 



(24) 

To calculate NTMEs M a of the (/3+/3+) 0l/ and (e/? + ) 0l , 
modes, the matrices [fz,N] a p an d [Fz.jv are eval- 

uated using expressions given by Eqs. (|2"3")l and 
respectively. The required NTMEs M a are obtained us- 
ing Eq. (|20p with 20 gaussian quadrature points in the 
range (0, 7r). 



III. RESULTS AND DISCUSSIONS 

The model space, single particle energies (SPE's) and 
parameters of the effective two-body interaction are the 
same as our earlier calculations on (e + /3/3) 2l/ decay of 



96 Ru, 102 Pd, 106,108 Cd [23^ 124,126 Xe5 130,132 B& g| &nd 

156 Dy [13] isotopes for the 0+ — > 0+ transition. We 
briefly present a discussion about them for the sake of 
completeness as well as present convenience. The dou- 
bly even 76 Sr (N = Z = 38) and 100 Sn (N = Z = 50) 
nuclei were treated as inert cores for the nuclei in the 
mass region A = 96 — 108 and A = 124 — 156, respec- 
tively. The change of model space was forced upon be- 
cause the number of neutrons increase to about 40 for 
nuclei occurring in the mass region A = 130 and with 
the increase in neutron number, the yrast energy spectra 
was compressed due to increase in the attractive part of 
effective two-body interaction. In Table U we have given 
the single particle orbits, which span the valence space 
and corresponding SPEs. In the model space with 76 Sr 
core, the lpi/2 orbit was included to examine the role of 
the Z — 40 proton core vis-a-vis the onset of deformation 
in the highly neutron rich isotopes. For 156 Dy and 156 Gd 
isotopes, the SPE's used for 0hu/ 2 , I/7/2 an d 0h 9 / 2 or- 
bits were 4.6 MeV, 11.0 MeV and 11.6 MeV, respectively. 



TABLE I: Single particle orbits of the model space and SPEs 
for protons and neutrons. 



A = 96 


- 108 


A = 124 


- 156 


Orbits 


e (MeV) 


Orbits 


e (MeV) 


lPl/2 


-0.8 


2si/ 2 


1.4 


2Sl/ 2 


6.4 


1^3/2 


2.0 


l<*3/2 


7.9 


l<*5/2 


0.0 


Ids/a 


5.4 


1/7/2 


12.0 


057/2 


8.4 


0fl7/2 


4.0 


0ff9/2 


0.0 


Ohg/2 


12.5 


0ftll/2 


8.6 


0hu/2 


6.5 



The HFB wave functions were generated by using an 
effective Hamiltonian with PQQ type of effective two- 
body interaction [33| given by 



H = H sp + V(P) + ( qq V(QQ), 



(25) 



where H sp , V(P) and V(QQ) represent the single parti- 
cle Hamiltonian, the pairing and quadrupole-quadrupole 
part of the effective two-body interaction, respectively. 
The arbitrary parameter (^ qq was introduced to study the 
role of deformation by varying the strength of QQ in- 
teraction and the final results were obtained by using 
Q q = 1. Following Heestand el al. [IH, who have used 
G p = 30 1 A MeV and G n = 20 /A MeV to explain the 
experimental <?(2 + ) data of some even-even Ge, Se, Mo, 
Ru, Pd, Cd and Te isotopes in Greiner's collective model 
[III, we used the same strengths for A = 96 — 108 nuclei. 
In the case of A = 124—132 isotopes, the strengths of the 
pairing interaction were fixed as G p = G n = 35/^4 MeV. 
However, we used G p = G n = 30 /A MeV for 156 Dy and 
156 Gd isotopes. 

The parameters of the QQ interaction were fixed as fol- 
lows. The strengths of the like particle components Xpp 
and Xnn were taken as 0.0105 MeV 6~ 4 , where b is os- 
cillator parameter. The strength of proton-neutron (pn) 
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component \pn was varied so as to obtain the spectra of 
considered nuclei A = 96 — 156 in optimum agreement 
with the experimental data. The theoretical spectra was 
taken to be the optimum one if the excitation energy of 
the 2 + state E 2 + was reproduced as closely as possi- 
ble to the experimental value. All the parameters were 
kept fixed throughout the subsequent calculations. The 
reliability of HFB wave functions was tested by obtaining 
an over all agreement between theoretically calculated re- 
sults for the yrast spectra, reduced B(E2:Q + — > 2 + ) tran- 
sition probabilities, static quadrupole moments Q{2 + ) as 
well as (/-factors g{2 + ) of the above mentioned nuclei and 
the available experimental data. The same PHFB wave 
functions were employed to calculate NTMEs M 2v and 
half-lives T^ 2 (0+ -> 0+) of {e + /3(3) 2v decay for 96 Ru, 

102 Pd, 106,108 Cd j2^ ) 124,126 Xej 130,132 Ba j^J and 156 Dy 

[H| isotopes. It was also shown that the proton-neutron 
part of the PQQ interaction, which is responsible for trig- 
gering deformation in the intrinsic ground state, plays an 
important role in the suppression of M% v . 



A. Results of (P + I3 + ) 0v and (e/3 + ) 0l , modes 

The phase space factors Gqi of (P + P + ) 0l/ and (e/3 + ) , y 
modes have been evaluated by Doi et al. with gA = 1.261 
@. We use the phase space factors after reevaluat- 
ing them for gA = 1.254. The phase space factors of 
f3 + (3 + (ef3 + ) modes (in yr _1 ) used in the present calcu- 
lation are 2.243 x 1(T 18 (2.664 x 1CT 17 ), 2.532 x 1CT 18 
(3.635 x 1CT 17 ), 3.048 x 1CT 18 (5.654 x 1CT 17 ) and 
5.114 x 1(T 19 (4.901 x 10~ 17 ) for 96 Ru, 106 Cd, 124 Xe and 
130 Ba nuclei, respectively p. For 102 Pd and 156 Dy nu- 
clei, we calculate Gqi following the notations of Doi et al. 
@ in the approximation C\ — 1.0, C 2 — 0.0, C3 = 0.0 
and Ri.i{e) = R+i(e) + R-i(e) =1.0. The calculated 
G01 of the e(3 + mode for 102 Pd and 156 Dy isotopes are 
6.0 x 10~ 19 yr -1 and 3.250xl0 -17 yr _1 , respectively. 

In Table HQ the NTMEs Mp, M GT , M Fh and 
Mark required to study the (P + (3 + ) 0l/ and (e/3 + ) „ 
modes of 96 Ru, 102 Pd, 106 Cd, 124 Xe, 130 Ba and 156 Dy nu- 
clei are compiled. Following Haxton's prescription [Ztl . 
the average energy denominator is taken as A — 1.2A 1 ' 2 
MeV. We calculate the four NTMEs in the approxima- 
tion of point nucleons, point nucleons plus Jastrow type 
of SRC with Miller and Spencer parametrization [4l|, fi- 
nite size of nucleons with dipole form factor and finite 
size plus SRC. In the case of point nucleons, the NTMEs 
Mp and Mqt are calculated for A and A/2 in the en- 
ergy denominator. It is observed that the NTMEs Mp 
and Mqt change by 7.8-9.8% for A/2 in comparison to 
A in the energy denominator. Therefore, the dependence 
of NTMEs on average excitation energy A is small and 
the closure approximation is quite good in the case of 
(/3 + (3 + ) aiJ and (e/3 + ) „ modes as expected. In the ap- 
proximation of light neutrinos, the NTMEs Mp and Mqt 
are reduced by 17.8-21.4% and 12.2-14.2% for point nu- 
cleoli plus SRC, and finite size of nucleons respectively. 



Finally, the NTMEs change by 21.7-25.8% with finite 
size plus SRC. In the case of heavy neutrinos, the Mph 
and M GTh get reduced by 33.9-38.0% and 65.0-68.5% 
with the inclusion of finite size and finite size plus SRC. 



0.5 n 




FIG. 1: Radial dependence of Cf{t), Cgt{t) and Co v (r) with 
FNS and SRC effects for the (f3 + f3 + ) Q and (e/3+) decay 
modes of 96 Ru, 102 Pd, 106 Cd, 124 Xe, f30 Ba and 156 Dy iso- 
topes. 
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TABLE II: Calculated NTMEs for the 0+ -> 0+ transition of {P + P + ) 0l/ and {sP + ) 0l/ modes in the mass mechanism. 



Nuclei 


NTMEs 


Point 

A 


A/2 


Point+SRC 


Extened 


Extended+SRC 


96 Ru 


M F 


0.4983 


0.5372 


0.3969 


0.4309 


0.3757 




1V1QT 


-2.4780 


-9 6896 


-9 0000 


-2.1591 


-1.8992 






35 891 7 




o 


22.4117 


11.4829 




iviQTh 


-1f!9 391 




Q 


-106 353 


-54.7130 


102p d 


M F 


0.6464 


0.6995 


0.5233 


0.5632 


0.4965 






-2.7663 


-2.9861 


-2.1863 


-2.3785 


-2.0631 




Mtpu 
lvl Fh 


43.3140 




o 


28.1494 


14.7508 




iviQTh 


-904 336 




q 


-129.721 


-67.0114 


106 Cd 


M F 


0.9583 


1.0394 


0.7704 


0.8319 


0.7299 




IVIQT 


-4.3495 


-4.7284 


-3.4635 


-3.7594 


-3.2769 






66 1 1 96 




o 


49 5Q8Q 


99 1 888 






-31 1 922 




o 


-197.061 


-101.408 


124 Xe 


M F 


0.4865 


0.5333 


0.3915 


0.4233 


0.3717 






-2.1387 


-2 3299 


-1.6905 


-1.8416 


-1.5978 




ilfp, 
±\±Fh 


33 7569 




o 


21.1455 


1 8449 

1W . O -L lis 




ivJ-GTh 


-159 250 




o 


-98 9817 


-50.4944 


130 Ba 


M F 


0.4183 


0.4593 


0.3338 


0.3623 


0.3163 




M gt 


-1.8626 


-2.0325 


-1.4633 


-1.5986 


-1.3812 




M Fh 


30.0461 







18.7025 


9.5438 




Moth 


-141.744 







-87.8418 


-44.6828 


156 Dy 


M P 


0.2461 


0.2698 


0.2022 


0.2160 


0.1926 




M GT 


-1.1281 


-1.2319 


-0.9208 


-0.9867 


-0.8754 




M Fh 


15.7014 







10.3729 


5.4997 




Moth 


-74.0722 







-48.6696 


-25.6980 



The radial dependence of Cq u {t) defined by 

oo 

M „= fc 0v (r)dr (26) 
o 

has been studied in the QRPA by Simkovic et al. [39[ 
and ISM by Menendez et al. In both QRPA and 

ISM calculations, it has been established that the con- 
tributions of decaying pairs coupled to J = and J > 
almost cancel beyond r w 3 fm and the magnitude of 
Coy(r) for all nuclei undergoing (f3~(3~) 0u decay are the 
maximum about the internucleon distance r ~ 1 fm. In 
Fig. [TJ we plot the radial dependence of the total matrix 
elements Cq v (r) as well as their Fermi and Gamow- Teller 
components due to the exchange of light neutrinos. It is 
noticed that the maximum value of CV(r), Cqt{t) and 
Coi/{r) is at r = 1.25 fm in agreement with the works 
done by Simkovic et al. [3^| and Menendez et al. [47| . 

In Table IIII1 we tabulate the extracted limits on the 
effective light neutrino mass < m„ > as well as heavy 
neutrino mass < Mn > using presently available exper- 
imentally observed limits on half- lives of (/3 + /3 + ) „ and 
(e/? + ) 0l/ modes. It is observed that limits on < m„ > 



and < Mn > are not so much stringent as in the case of 
(/3~/3~) 0l , mode. Further, better limits are obtained in 
the case of (£/3 + ) „ mode even for equal limits on half- 
lives of (/? + /3 + ) „ and (e/3 + ) 0l/ modes. In the case of 
{e(3 + ) 0l/ mode, the best limits obtained for 130 Ba nuclei 
are < m v >< 6.8 x 10 2 eV and < M N > > 2.25 x 10 4 
GeV. 

In Table IIV1 we compile available theoretical results 
in other nuclear models along with ours. To the best 
of our knowledge, no theoretical result and experimental 
half-life limit is available for 102 Pd and 156 Dy isotopes. 
Staudt et al. [52] have reported only NTMEs \M 0u \ = 
\Mqt — Mp\ in the mass mechanism. In the QRPA cal- 
culations of Hirsch et al. |l3| and Staudt et al. [52| . 
the former used two major oscillator shells, where as 
the latter used a model space consisting of 3fku + Atko + 
O/19/2 + 0hii/ 2 orbits. The used SPEs are identical. Both 
the calculation use a realistic effective two body interac- 
tion using Paris potential. The NTMEs \M$ V \ are almost 
identical in both the QRPA calculations but for 124 Xe, 
where a difference by a factor of 1.8 approximately is 
noticed. In the SQRPA model, Stoica et al. [53| have 
studied (/3 + /3 + ) „ and {e(3 + ) 0u modes of 106 Cd isotope 
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TABLE III: Upper and lower bounds on light and heavy neutrino effective masses < m„ > and < Mjv >, respectively, for the 
(/3+/3 + ) „ and (e/3 + ) 0l , modes of 96 Ru, 106 Cd, 124 Xe and 130 Ba isotopes. 



Nuclei 




^(yr) 

8(3+ 


Ref. 


< m v 


> (eV) 

sP+ 


< M N 
/3+/3+ 


> {GeV) 

Ef3+ 


yy Ru 


> 3.1 x 10 lb 


> 6.7 x 10 lb 


[48] 


8.52 x 10 b 


1.68 x 10 b 


16.38 


82.98 


106 Cd 


> 1.4 x 10 19 


> 7.0 x 10 19 


[49] 


2.14 x 10 4 


2.53 x 10 3 


6.90 x 10 2 


5.85 x 10 3 


124 Xe 


> 4.2 x 10 17 


> 1.2 x 10 18 


[50] 


2.29 x 10 5 


3.15 x 10 4 


65.12 


4.74 x 10 2 


130 Ba 


> 4.0 x 10 21 


> 4.0 x 10 21 


[51] 


6.66 x 10 3 


6.80 x 10 2 


2.30 x 10 3 


2.25 x 10 4 



TABLE IV: Predicted half-lives T"J 2 < m v > 2 °f {P + P + )q,j an d ( £ <^ + )oi/ mo< A es due to the exchange of light neutrino and 



extracted limits on effective heavy neutrino mass (Mn) 
denote WS and AWS basis respectively in reference [HEj . 



from the same predicted half- lives for < m„ >= 1 eV. The t and $ 



Nuclei Model 



Ref. 



Mr, 



Mr, 



\M a 



T[y 2 < m v 
(yr eV 2 ) 



Mr, 



IMn 



e/3" 1 



< Mjv > 
(GeV) 



Ru 



Pd 



Cd 



Xe 



J Ba 



PHFB 




0.376 


-1.899 


2.275 


MCM 


54 


-0.705 


1.678 


2.383 


QRPA 


13 


-0.98 


2.62 


3.60 


QRPA 


52 






4.228 






0.497 


-2.063 


2.560 


PHFB 




0.730 


-3.277 


4.007 


MCM 


54 


-1.191 


2.203 


3.394 


SQRPA(l) 


53 


-2.12 


5.73 


7.85 


SQRPA(s) 


53 


-2.18 


5.99 


8.17 


QRPA 


13 


-1.22 


3.34 


4.56 


QRPA 


52 






4.778 


PHFB 




0.372 


-1.598 


1.970 


MCM 


54 


-2.572 


5.729 


8.301 


QRPA 1 " 


55 


-2.236 


5.128 


7.364 


QRPA* 


55 


-2.574 


5.733 


8.307 


QRPA 


13 


-1.35 


3.92 


5.27 


QRPA 


52 






2.975 


PHFB 




0.316 


-1.381 


1.697 


MCM 


54 


-1.748 


3.382 


5.130 


QRPA 


13 


-1.50 


4.02 


5.52 


QRPA 


52 






5.579 



2.249xl0 28 
2.050xl0 28 
8.981xl0 27 
6.511xl0 27 



6.424xl0 27 
8.953xl0 27 
1.674xl0 27 
1.545xl0 27 
4. 960 xlO 27 
4.517xl0 27 

2.208xl0 28 
1.243xl0 27 
1.580xl0 27 
1.241xl0 27 
3.084xl0 27 
9.678xl0 27 



1.772xl0 2 
1.940xl0 2 
1.676xl0 2 
1.641xl0 2 



1.894xl0 27 
1.726xl0 27 
7.563 xlO 26 
5.483xl0 26 

6.643 xlO 28 



474 xlO 26 
236 xlO 26 
166 xlO 26 
076 xlO 26 
455 xlO 26 
146 xlO 26 



1. 191xl0 2 ' 
6. 703 xlO 25 
8. 517xl0 25 
6. 693 xlO 25 
1.663 xlO 26 
5. 218 xlO 26 

1.849xl0 27 
2.025xl0 26 
1.749xl0 26 
1.712xl0 26 



11.483 -54.713 66.196 1.40xl0 7 



14.751 -67.011 81.762 1.53xl0 7 



22.189 -101.408 123.597 1.48xl0 7 



10.845 -50.494 61.339 



9.544 



1.49x10' 



-44.683 54.227 1.53xl0 7 



3 Dy 



PHFB 



0.193 -0.875 



1.068 



7.044xl0 2 



5.500 



-25.698 31.198 



1.40x10' 



using two model spaces, namely small basis (oscillator 
shells of 3Huj — 5Uuj + «i 3 / 2 orbits) and a large basis (os- 
cillator shells of 2huj — 5Hlo + H3/2 orbits) with two-body 
effective interactions derived from the Bonn-A potential. 
The NTMEs calculated in the SQRPA [H do not de- 
pend much on the model space and differ by a factor of 
1.8 approximately from those of Hirsch et al. [13j. In the 
MCM, Suhoncn et al. ^ have studied the {/3 + /3 + ) 0v 
and (£/?+) „ modes of 96 Ru, 106 Cd, 124 Xe and 130 Ba nu- 
clei. It is worth mentioning that besides the model space, 
SPEs and effective two-body interaction, different values 
of g A , specifically g A = 1.254 [HiiHI and L0 HIHl: 



are also used in these calulations. 

The calculated NTMEs \M 0v \ in the PHFB model for 
the 96 Ru and 106 Cd isotopes are very close to those ob- 
tanied in the MCM, and in the later case also to the 
QRPA results. For 124 Xe and 130 Ba isotopes, the NTMEs 
are smaller than those in other models and this is re- 
flected in half-lives which are up to one order of mag- 
nitude longer. As the extracted limits on the effective 
neutrino masses < m v > and < Mn > are not strin- 
gent enough, it is more meaningful to calculate half-lives 
of (/3 + /3 + ) „ and (e/3 + ) 0l , modes, which will be useful 
for the design of future experimental set ups. Hence, 
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we calculate half-lives of ((3 + (3 + ) Qu and {e(3 + ) 0v modes 
for < m v >= 1 eV and extract corresponding limits on 
heavy neutrino mass < Mn >, which are given in the 
same Table [TV] 

In the mass mechanisms, there are two noteworthy ob- 
servations. The equality in NTMEs of {(3 + (3+) av and 
(e/? + ) 0[/ modes implies that 



T°; 2 (/?+/?+)_ G 01 (£/?+) 



T^ 2 (e{3+) G 01 ((3+(3+Y 



(27) 



Therefore, the experimental observation of (e/3 + ) 0l / mode 
will provide the half-life T°fo {(3+ (3+) of {(3 + (3+) Q „ mode 
as the phase space factors are exactly calculable. Fur- 
ther, it is noticed that the ratios of \Mq„\ and |Moat| 
given in Table |TT] are almost constant for different nu- 
clei and | Moat | / |Mo„| « 29 — 32 approximately. Similar 
behaviour of the ratios | AfoAf | / \Mq v \ ~ 28 — 30 is also 
observed for the NTMEs of ((3~(3~) 0u mode [25]. This 
implies that in the mass mechanism, the half-lives for dif- 
ferent nuclei due to exchange of light and heavy neutrinos 
are also in constant ratio 



IM, 



0N\ 



\M 0v \ 



(28) 



It will be interesting to verify whether the observed con- 
stancy of |Mojv| / \Mq v \ in different nuclei is a generic 
feature or artifact of the present calculation. 



B. Quadrupolar correlations and deformation 
effects 



As already mentioned, the quadrupolar correlations 
are mainly responsible for the deformation of nuclei. 
To understand the role of deformation on NTMEs M a 
{a = F, GT, Fh, GTh) of {[3+ (3 + ) 0u and (e/3 + ) 0l/ modes, 
we investigate the variation of the latter by changing the 
strength of the QQ interaction ( qq for the case in which 
NTMEs are calculated with finite size and short range 
correlations. It is observed that in general, there is an 
inverse correlation between the magnitudes of NTMEs 
and quadrupole moments Q(2 + ) as well as deformation 
parameters /?2- Further, the effect of deformation on M a 
is quantified by defining a quantity D a as the ratio of 
M a at zero deformation ((^ qq — 0) and full deformation 
(( qq = 1). The D a is given by 



M a (C qq = 0) 

M a (( qg = iy 



(29) 



The tabulated values of D a in Table |V] for 96 Ru, 102 Pd, 
106 Cd, 124 Xe, 130 Ba and 156 Dy nuclei suggest that the 
NTMEs M a are suppressed by factor of 1.7-10.7 in the 
mass range A = 96 — 156 due to deformation effects. We 
also give the same deformation ratio Di v for comparison 
in the last row of the same table, which also change by al- 
most same amount due to the deformation effects. Hence, 



it is clear that the deformation effects are important for 
((3 + (3 + ) Qv and (e/3 + ) „ modes as well as (e + (3f3) 2u decay 
so far as the nuclear structure aspect of e + (3(3 decay is 
concerned. 

In the left and right panels of Fig. [2] and [3j we present 
the variation of NTMEs \M Qv \ and \M 0N \ due to the light 
and heavy neutrino exchange, respectively, with respect 
to = /^(parent)— ^{daughter) for the above men- 
tioned e + (3(3 emitters. The theoretically calculated de- 
formation parameters for parent and daughter nuclei 
have been given in Refs. [23|, [24[ and we present them 
in Table IVII for convenience. It can be noticed that the 
variation in |Moj,| with changing A/3a is similar as that of 
I Moat I . Moreover, it can be observed in Fig. [2] and [3] that 
the NTMEs remain constant even when one of the nuclei 
is spherical or slightly deformed. With further increase 
in deformation, the NTMEs in general become the maxi- 
mum for A/?2 = and then decrease with increase in the 
difference between the deformation parameters. To sum- 
marize, the independent deformations of initial and final 
nuclei are important parameters to describe the NTMEs 
Mq v and Mqat of {(3 + (3 + ) 0u and (£(3 + ) 0u modes. 



TABLE V: Ratios D a for 96 Ru, 102 Pd, 106 Cd, 124 Xe, 130 Ba 
and 156 Dy isotopes. 



Ratios 


*Ru 




iub Cd 


m Xe 


13U Ba 


lbb Dy 


D F 


2.92 


2.52 


1.91 


3.83 


4.68 


10.42 


Dgt 


2.48 


2.73 


1.96 


3.88 


4.72 


10.68 


D Fh 


2.61 


2.34 


1.72 


3.42 


4.11 


10.20 


DGTh 


2.49 


2.36 


1.72 


3.45 


4.13 


10.20 


D 2v 


3.13 


3.40 


2.06 


3.63 


4.66 


13.64 



TABLE VI: Calculated [H, HI 
mation parameters (32 of parent 
pating in {(3 + (3 + ) Qv and (e(3+) 



and experimental [5g| defor- 
and daughter nuclei partici- 
modes. 



Nuclei 


Theory 


(32 

Experiment 


y6 Ru 


0.161 


0.1579±0.0031 


96 Mo 


0.191 


0.1720±0.0016 


102 pd 


0.185 


0.196±0.006 


102 Ru 


0.232 


0.2404±0.0019 


106 Cd 


0.176 


0.1732±0.0042 


106 Pd 


0.203 


0.229±0.006 


124 Xe 


0.210 


0.212±0.007 


i24 Te 


0.164 


0.1695±0.0009 


130 Ba 


0.234 


0.2183±0.0015 


130 Xe 


0.166 


0.169±0.007 


156 Dy 


0.300 


0.2929±0.0016 


156 Gd 


0.316 


0.3378±0.0018 
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% Ru ^% Mo % Ru^ % Mo 




FIG. 2: NTMEs of (/?+ P + ) and (e/3 + ) Qv modes for 96 Ru, 102 Pd, 106 Cd isotopes due to the exchange of light (left hand side) 
and heavy (right hand side) neutrinos as a function of the difference in the deformation parameters A/32, "x" denotes the 
NTME for calculated A/3 2 at Q qq = 1. 



IV. CONCLUSIONS 

We have calculated the NTMEs M F , M GT , M Fh and 
Mark required to study the {P + (3 + ) 0v mode of 96 Ru, 
106 Cd, 124 Xe and 130 Ba as well as the (e(3+) 0lJ mode 
of 96 Ru, 102 Pd, 106 Cd, 124 Xe, 130 Ba and 156 Dy nuclei 
for the + — * + transition in the Majorana neutrino 
mass mechanism using the set of HFB wave functions, 
the reliability of which was tested by obtaining an over- 
all agreement between theoretically calculated results 
for the yrast spectra, reduced B(E2:0 + — > 2 + ) transi- 
tion probabilities, static quadrupole moments Q(2 + ) and 
g-factors g(2 + ) and NTMEs Mi v as well as half- lives 
T^J 2 of (e + /3/3) 2i , decay and the available experimental 
data [H, [U |34|. The existing experimental data on 
(/3 + /3 + ) „ and (e/3 + ) 0l/ modes fail to provide stringent 
limits on the extracted effective mass of light neutrino 



(m v ) and heavy neutrino (Mjv). Hence, we calculate 
half-lives Ty 2 of these modes for the light neutrino and 
extract limits on (Mn). In the mass mechanism, the half- 
lives Ty 2 (/3+P+) and T°J 2 (e/3+) are related through the 
exactly calculable phase space factors Goi (/3 + /3 + ) and 
Gqi In addition, it is observed that the ratio 

of NTMEs | M 0A r | / \M 0v \ w 30 is a constant for differ- 
ent nuclei so that half-lives due to the exchange of light 
and heavy neutrinos are also in constant ratio. Further, 
the role of deformation on NTMEs M F , M GT , M Fh and 
Mark for (/3 + /3 + ) „ and (£f3 + ) 0l/ modes is investigated 
by changing the strength Q qq of the QQ interaction. It is 
noticed that there is an inverse correlation between the 
magnitudes of NTMEs and quadrupole moments Q(2 + ) 
as well as deformation parameters 02 ■ The NTMEs are 
suppressed by factors of 1.7-10.7 in the considered mass 
range A — 96 — 156 implying that the nuclear struc- 
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ture effects are also important for (P + (3 + ) 0l/ and (e/3 + ) 0l , 
modes. The deformation of individual nucleus is an im- 
portant parameter for calculating NTMEs Mo„ and Moat 
of (/3 + /? + ) „ and (e/3 + ) 0iy modes. 
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